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ABSTRACT 

In this paper we report the results of axisymmetric relativistic MHD simulations 
for the problem of Kerr black hole immersed into a rarefied plasma with "uniform" 
magnetic field. The long term solution shows properties which are significantly dif- 
ferent from those of the initial transient phase studied recently by Koide(2003). The 
topology of magnetic field lines within the ergosphere is similar to that of the split- 
monopole model with a strong current sheet in the equatorial plane. Closer inspection 
reveals a system of isolated magnetic islands inside the sheet and ongoing magnetic 
reconnection. No regions of negative hydrodynamic " energy at infinity" are seen inside 
the ergosphere and the so-called MHD Penrose process does not operate. Yet, the ro- 
tational energy of the black hole continues to be extracted via purely electromagnetic 
mechanism of Blandford and Znajek(1977). However, this is not followed by devel- 
opment of strong relativistic outflows from the black hole. Combined with results of 
other recent simulations this signals a potential problem for the standard MHD model 
of relativistic astrophysical jets should they still be observed at distances as small as 
few tens of gravitational radii from the central black hole. 

Key words: black hole physics - magnetic fields - methods:numerical. 



1 INTRODUCTION 

Observations of various astrophysical phenomena such as 
active galactic nuclei and galactic microquasars often reveal 
powerful relativistic jets streaming away from a massive cen- 
tral object, most likely a black hole. The exact mechanism of 
generating powerful relativistic jets from black holes system 
is not yet known, although a number of interesting models 
have been put forward and studied with various degree of 
detail during the last few decades. 

By now there has emerged a general consensus that the 
central engine should involve a rotating black hole linked 
with the jets by means of strong magnetic field and an ac- 
cretion disc supporting this magnetic field by its electric 
currents. This magnetic field is believed to serve a number 
of key functions: 1) to power the jets via extracting the rota- 
tional energy of the black hole; 2) to provide the jet collima- 
tion via magnetic hoop stress; and 3) to suppress mixing of 
the rarefied jet plasma with the relatively dense surround- 
ing medium such as the coronal plasma of the accretion disc, 
which is needed to ensure that the jet plasma remains mag- 
netically dominated and, thus, can be accelerated to the 
ultra-relativistic speeds inferred from the observations. 



The most simplified mathematical frameworks for 
highly magnetised relativistic plasma are an ideal relativis- 
tic magnetohydrodynamics (MHD) and magnetodynamics 
(MD, which can be described as MHD in the limit of zero 
particle inertia, e.g Komissarov 2002). Both frameworks, 
however, involve rather complex equations particularly in 
the case of curved spacetime. This explains a relatively slow 
progress in the theory - only a limited number of analytical 
solutions have been found so far and only for the limiting 
case of a slowly rotating black hole where one can employ 
the perturbation method. The most important of them is the 
MD solution for a monopole magnetosphere by Blandford 
and Znajek (1977). Indeed, this solution describes an outgo- 
ing Poynting flux from a Kerr black hole as if the hole was 
a magnetized rotating conductor whose rotational energy 
can be efficiently extracted by means of magnetic torque. 
This most important feature of the Blandford-Znajek solu- 
tion (BZ) is also the most intriguing and puzzling as such 
conductor does not really exist. 

Often, this missing element of the BZ model is ar- 
tificially introduced in the form of the so-called "mem- 
brane", or the "stretched horizon", located somewhat 
above the real event horizon (Macdonald & Thome 1982 
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IThorne et al.l986ll . This "physical" interpretation of the BZ 
solution provides simple means of communicating the re- 
sults to wide astrophysical community and for this reason it 
has been almost universally accepted. However, its artificial 
nature and the fact that no proper physical interpretation 
had been pushed forward also made the result vulnerable to 
criticism on theoretical grounds and stimulated attempts to 
find alternative ways of magnetic extraction of rotational en- 
ergy of black holes jPunsly fc Coroniti 1990a] |Punsly 200l| 
Takahashi ct al.1990}. These attempts seem to draw the in- 
spiration from a completely non-magnetic mechanism pro- 
posed much earlier by Penrose (1969). The key role in this 
mechanism is played by the ergosphere of a rotating black 
hole. Within the ergosphere particles can acquire negative 
energy (or rather "energy at infinity") and if such particles 
are swallowed by the black hole its energy can decrease. In 
the original Penrose process such particles are created via 
close range interaction (collisions, decay) with other parti- 
cles which gain positive energy and carry it away. However, 
electrically charged particles can also be pushed onto or- 
bits with negative energy by the Lorentz force and this is 
what makes possible the so-called "MHD Penrose process" 
IjTakahashi et al. 19901 IKoide et al.2002l IKoide 20031 which 
is the common element of all alternative magnetic mecha- 
nisms. 

Recently there has been a renewed interest to 
this problem. The main reason for this is the ar- 
rival of robust numerical methods for relativistic 
astrophysics, e.g. (IPons et al. 19981 IKomissarov 1999 
IKoide et al.19991 IKoldoba et al.27)?)2l iGammie et al.2003l 
IDel Zanna et al.2003l |De Villiers fc Hawley 20051 . Time- 
dependent MD and MHD simulations of monopole magne- 
tospheres of black holes have demonstrated the asymptotic 
stability of the BZ solution (Komissarov 2001 1 and provided 
additional arguments in favour of the MD approximation in 
this particular case (Komissarov 2004a Komissarov 2004b). 
Moreover, simulations of magnetically driven accretion 
disks have shown the development of low density axial 
regions, "funnels", closely described by the BZ solution 
l|De Villiers et al. 20031 |McKinney fc Gammie 2004) . 

On the other hand, the MHD simulations of a black 
hole immersed into a rarefied plasma with uniform magnetic 
field seemed to provide support for the MHD Penrose model 
IIKoide et al.2002l IKoide 20031 . At least, it was found that 
by the end of the simulations the inflow of particles with 
negative energy at infinity accounted for about one half of 
the extracted energy. Unfortunately, the termination time 
of these simulations was surprisingly short, approximately 
one half of the rotational period of the black hole (presum- 
ably due to some computational problems.) As the result 
one could not tell whether the MHD Penrose mechanism 
would remain effective on a long-term basis. In this paper 
we present the results of new simulations which were carried 
out for much longer period of time and which bring new light 
on this important issue. 



2 BASIC EQUATIONS AND NUMERICAL 
METHOD 

In these simulations we solve numerically the equations of 
ideal MHD in the space-time of a Kerr black hole. This 



space-time is described using the foliation approach where 
the time coordinate t parametrises a suitable filling of 
spacetime with space-like hypersurfaces described by the 3- 
dimensional metric tensor 7^. These hypersurfaces may be 
regarded as the "absolute space" at different instances of 
time t. If {x 1 } are the spatial coordinates of the absolute 
space then the metric form can be written as 



ds = (/? — ol )dt + 2f3idx l dt + 'yijdx 1 dx J 



(1) 



where a > is called the "lapse function" and (3 is the 
"shift vector" . In this study we employ the Kerr-Schild coor- 
dinates, t, (f>, r, 9, which like the well known Boyer-Lindquist 
coordinates ensure that none of the components of the met- 
ric form depend on t and 0. At spatial infinity these co- 
ordinate systems do not differ but the Kerr-Schild system 
does not have a coordinate singularity on the event hori- 
zon. Other details can be found in e.g. (Komissarov 2004a 
|McKinney fc Gammie 2004) . 

The evolution equations of ideal MHD include the con- 
tinuity equation, 



d t (ct^jpu) + d l (a y ^ypu) = 0, 
the energy-momentum equations, 

and the induction equation, 
8t(B t )+e^ k d J (E k )=0. 



(2) 



(3) 



(4) 



Here g a f3 is the metric tensor of spacetime, 7 = det^ij), 
e %ik is the Levi-Civita pseudo-tensor of space, p is the proper 
mass density of plasma, u" is its four-velocity vector. The 
total stress-energy-momentum tensor, is a sum of the 

stress-energy momentum tensor of matter, 



r > = wu u 

(m) 



pg 



(5) 



where p is the thermodynamic pressure and w is the enthalpy 
per unit volume, and the stress-energy momentum tensor of 
electromagnetic field, 



(6) 



where F vtl is the Maxwell tensor of the electromagnetic field. 
The electric field, E, and the magnetic field, B, are defined 
via 



Ek = 
and 



— p ■ 1 

2 l 3 k ' 



B % — a *F l 



(7) 



(8) 



where "F^ is the Faraday tensor of the electromagnetic 
field, which is simply dual to the Maxwell tensor. In the 
limit of ideal MHD 



E = vxB or Ei — eijuv 3 B 



3 n k 



(9) 



where v l — u 1 /u is the usual 3- velocity of plasma. In ad- 
dition, the magnetic field must satisfy the divergence free 
condition 



0. 



(10) 



Note, that 1) all the components of vectors and tensors ap- 
pearing in these equations are the those measured in the 
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Figure 1. The large scale structure of the numerical solution 
during the initial phase, t = 6. The colour image shows the dis- 
tribution of rest mass density ,log 10 p. The contours show the 
magnetic flux surfaces (the lines of the poloidal component of 
magnetic field.) The arrows show the poloidal component of ve- 
locity relative to the coordinate grid. 

coordinate basis, {&,}, of the Kerr-Schild coordinates; 2) 
throughout the paper we employ such units that the speed 
of light c = 1, the gravitational constant G = 1, the black 
hole mass M = 1, and the factor <in does not appear in 
Maxwell's equations. 

Our numerical scheme is a 2D Godunov-type upwind 
scheme which utilises special relativistic Riemann solver de- 
scribed in Komissarov(1999) and uses the method of con- 
straint transport <Evans fc Hawley 1988) to preserve the 
magnetic field divergence free. Other details are outlined 
in lIKomissarov 2 004b I . 



3 NUMERICAL SIMULATIONS 
3.1 Setup 

In these simulations the rotational parameter a of the Kerr 
metric is a = 0.9 which gives the event horizon radius r+ ~ 
1.44. The axisymmetric computational domain covers < 
6 < n and 1.35 < r < 53.1. The computational grid has 401 
cells in the redirection, where it is uniform, A9 — const , and 
400 cells in the r-direction. The cell size in the r-direction, 
Ar, is such that the corresponding physical lengths in both 
directions are equal in the equatorial plane. 

The usual axisymmetric boundary conditions are im- 
posed at 9 = and 6 — n boundaries. At the outer bound- 
ary, r — 53.1, the initial values of all variables are imposed 
throughout the whole run, the termination time for the sim- 
ulations was set to t — 60 so no waves emitted from the 
dynamically active region near the black hole had a chance 



to get reflected of the boundary and effect the inner solution. 
The inner boundary r = 1.35 is well inside the event horizon, 
which justifies use of "radiative boundary conditions" . 

The initial velocity field is set to be the same as the one 
of the fiducial observers of the Kerr-Schild foliation, who 
spiral towards the black hole (e.g. Komissarov 2004a). 

The initial electromagnetic field has the same "uniform" 
magnetic component, aligned with the rotational axis of the 
black hole, as in the vacuum solution of Wald: 

Fjiu = [m^u] + 2afc [MiI/ ]), (11) 

where k v — d t and m v = 8$ are the Killing vectors of the 
Kerr spacetime llWald 1974t . However, the initial electric 
component is different as it has to satisfy the condition of 
perfect conductivity 

The initial thermodynamic pressure and the rest mass 
density of plasma are set to be p = 0.01p m and p — 0.05p m , 
where p m = {B 2 — E 2 )/2 is the magnetic pressure. Finally, 
the equation of state employed in the simulations describes 
polytropic gas with the ratio of specific heats, F = 4/3. 

Summarising, although the setup of our simulations is 
not exactly the same as in Koidef2002} it is quite similar. 
In both cases we are dealing with magnetically dominated 
plasma of similar magnetisation. In both cases the initial 
magnetic field is described by the Wald vacuum solution 
( Wald 1974). Thus one would expect to obtain at least qual- 
itatively similar solution in the common region of computa- 
tional domains. Because of different space-time foliations, 
Koide(2003) used the Boyer-Lindquist coordinates, the dif- 
ference in computational domains is not simply reduced to 
the difference in the range of the radial coordinate, r. Differ- 
ent definitions of the global time coordinate, t, should also 
be taken into account. This, however, is only significant in 
the spacial regions very close to the event horizon where the 
Boyer-Lindquist system becomes singular. 

All conservative schemes for the relativistic magnetohy- 
drodynamics have an upper limit on plasma magnetization 
above which they fail. At this limit, which somewhat varies 
from problem to problem and also depends on the resolu- 
tion, the numerical error for the total energy density be- 
comes comparable with the energy density of matter. This 
forces us to pump in fresh plasma in regions where its mag- 
netization becomes dangerously high. Since in such regions 
the dynamical role of particles is rather insignificant this 
measure should not have a strong effect in most respects. 
The critical condition we set in these simulations is 

wW 2 - p = 0.0FB 2 , (12) 

where W is the Lorentz factor of the flow and B is the mag- 
netic field strength as measured by the local FIDO. Should 
the energy density of matter drop below 0.01-B 2 , both p and 
p are artificially increased by the same factor. To minimise 
the effect of the mass injection on the flow the velocity of 
the injected matter is set to be equal to the local velocity 
of the flow. In fact, new particles must be constantly cre- 
ated in real magnetospheres of black holes but the details 
of this process can be rather different JBeskin et al.l99T1 
IHirotani fc Okamoto 19981 |Phinney 1982) . An additional 
lower limit was set on the value of the thermodynamic pres- 
sure, which was not allowed to drop below O.Olp. 
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Figure 2. Left panel: The ut component of plasma four-velocity and magnetic flux surfaces at t = 6. Right panel: the same as in the 
left panel but at t = 60. In both panels there is shown the same set of magnetic flux surfaces. The thick solid line shows the ergosphere. 
Note that —ut gives the specific "energy at infinity" for a free particle and that the region of negative hydrodynamic energy at infinity 
is even somewhat narrower than the region of positive ut because of the pressure contribution. 



3.2 Results and Discussion 

We have found that numerical solution exhibits two phases 
with rather different properties: 1) the rather short initial 
phase which is dominated by a rapid evolution in the neigh- 
bourhood of the black hole and 2) the final phase where 
solution settles to an approximate steady state in this re- 
gion. 

In the initial phase our numerical solution is indeed 
very similar to the one described in Koide (2003.3 • Plasma 
mainly slides along the magnetic field lines towards the 
equatorial plane where it passes through an accretion shock 
and forms an equatorial disc (see figO- The magnetic field 
lines are also pulled towards the black hole. Further away 
from the hole this seems to be mainly due the velocity 
field of the initial solution as the subsequent increase of 
the magnetic pressure in the central region quickly halts 
this motion and the magnetic field lines of these distant re- 
gions eventually straighten up (see fig0J. However, near the 
horizon this pulling is certainly caused by the hole as the 
field lines never straighten up there. Since similar effect is 
also observed in simulations of inertia-free magnetospheres 
(Komissarov 2004a) it cannot not be explained simply by 
dragging of the magnetic field lines along with the accreted 
plasma and has to have a more general cause. In any case, 
these findings are in strong contrast with the exclusion of 
magnetic flux by rotating black holes found in axisymmetric 
vacuum solutions l)Wald 19741 IBicak fc Janis 1985^ . Those 
vacuum solutions were used in the past to argue low effi- 
ciency of the BZ mechanism in the case of rapidly rotating 



black holes. Our results shows that this argument is incorrect 
and vacuum solutions have to be used with more caution. 

The left panel of figure |2] shows the inner region of our 
solution at t — 6 in more detail. The colour image describes 
the distribution of the covariant time component of plasma's 
four- velocity, Ut, whereas the solid lines show the magnetic 
flux surfaces (the lines of poloidal magnetic field). For a 
free particle —Ut gives its specific energy at infinity. Since 
the hydrodynamic part of energy at infinity has the volume 
density 



aT\. 



^ = a(— wutu* +i 



(13) 



it can only be negative if ut > 0; because of the pressure 
contribution the region of negative hydrodynamic energy at 
infinity is somewhat narrower than the region of positive Ut ■ 
What is even more important is that the flux of hydrody- 
namic energy at infinity 



aT\ = —awutu 



(14) 



has the opposite direction to the flow velocity v t — u z /tt* , 
as required in the MHD Penrose process, if and only if u t is 
positive. Figure |21 shows that well in agreement with find- 
ings of tKoide et al. 20021 IKoi de 2003 ) the inner part of the 
ergospheric disc does indeed have positive Ut during the ini- 
tial phase. Moreover, the magnetic field has a very similar 
topology too - all magnetic field lines threading the ergo- 
spheric disc have a turning point in the equatorial plane and 
do not cross the event horizon. 

The reason for developing positive ut looks attractively 
simple. Due to the inertial frames dragging, all plasma enter- 
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Figure 4. The large-scale structure of the solution at t = 60. Top left panel: the rest mass density, p\ Top right panel: the gas temperature, 
log 10 (P/p); Bottom left panel: H t f,=(the total electric current through the circular loop r =const)/27r; Bottom left panel: the angular 
velocity of magnetic field lines Q/Q^, where is the angular velocity of the black hole. The solid lines show the magnetic flux surfaces 
and the arrows show the poloidal velocity relative to the coordinate grid of the Kerr-Schild coordinates. 



ing the black hole ergosphere is forced to rotate in the same 
sense as the black hole (to be more precise this is what is 
observed by a distant observer. In Kerr-Schild coordinates 
d<j>/dt may have both signs.) As the result a differential rota- 
tion inevitably develops along the magnetic field lines pen- 
etrating the ergosphere which causes twisting of these lines. 
The Alfven waves generated in this manner propagate away 



from the ergosphere to infinity and establish an outflow of 
energy in the form of Poynting flux. Because of the energy 
conservation the ergospheric plasma reacts by moving onto 
orbits of lower energy. Provided this plasma remains in the 
space between the event horizon and the ergosphere for long 
enough it may indeed end up having negative energy. 

In the problem under consideration the strong magnetic 
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Figure 3. The hydro-dynamic (dashed line) and the electromag- 
netic (solid line) components of the energy flux (— T r t ) through 
the event horizon at t = 60. 



field keeps the plasma of the ergospheric disc from falling 
into the black hole and allows it to gain negative energy (see 
figHJ- However, such configuration cannot be sustain forever 
within ideal MHD. Since, the energy is constantly extracted 
along the field lines of the equatorial disc the energy of the 
disc is constantly going down and no steady-state can be 
reached. At some point the magnetic configuration would 
have to change so that all magnetic field lines entering the 
ergosphere also penetrate the event horizon (Another possi- 
ble option could be a non-steady behaviour with magnetic 
field lines pulled in and out of the event horizon all the time.) 
As one can see in the right panel of figure|2] where we present 
the solution at t — 60, this is more or less what is observed 
in our simulations. At around t = 20 the ergospheric disc is 
fully swallowed by the hole and a strong current sheet de- 
velops in its place. As the result, the structure of magnetic 
field becomes similar to that of the split-monopole model 
(Blandford & Znajc k~l977| . This configuration persists till 
the very end of the simulations when the whole solution in 
the inner part of the spacial domain seems to reach a state 
of an approximate equilibrium. 

In fact, one can see in fig|5|at least two magnetic islands 
in the equatorial current sheet. Such structure is known be 
quite typical for the tearing mode instability (e.g. Priest & 
Forbes, 2000) and was suggested in the context of black 
hole magnetospheres in Beskin (2003). Here, the islands 
are formed during slow reconnection events resulting in the 
gradual escape of magnetic flux from the event horizon, the 
inner islands disappear into the hole whereas the outer ones 
move away from it and supply hot plasma for a thin out- 
flow sheath clearly seen in the temperature plot in fig2] 
This evolution is indeed very slow and might be related to 



a somewhat excessive capture of magnetic field line during 
the early stages. 

The most important feature of the solution after the 
restructuring of magnetic field is the total disappearance of 
the regions with positive ut (see the right panel of fig[2)- 
Thus, the MHD Penrose process does no longer operate 
in the black hole ergosphere. This is confirmed in figure 
which shows that the hydrodynamic flux of energy at infin- 
ity through the event horizon, -QT r t(m) , is everywhere neg- 
ative. However, the electromagnetic flux of energy at infin- 
ity through the event horizon is positive almost everywhere 
with the exception of a very thin equatorial belt where it can 
be slightly negative (figEJ • Thus, the pure electromagnetic 
Blandford-Znajek mechanism continues to operate. More- 
over, the total flux of energy at infinity through the event 
horizon is positive and, in spite of the fact that the plasma 
magnetization is many orders of magnitude lower than that 
expected in typical astrophysical conditions, the BZ mech- 
anism allows to extract the rotational energy of the black 
hole. 

One important feature of the initial phase which has 
been found in the previous simulations fKoide 2003} and 
which persists in the longer run is the lack of any notice- 
able plasma outflow from the black hole. The only excep- 
tion is a thin sheath of radius r ~ 3, most clearly seen in 
the top right panel of fig|H) which is supplied with relativis- 
tically hot plasma (p/p — 10 2 ) via the reconnection process 
in the ergospheric current sheet. The outflow in the sheath 
is rather slow with poloidal speed relative to the local FIDO 
never exceeding v p ~ 0.7c; it is most likely driven by the 
gas pressure. This finding is in striking contrast with the 
results of MHD simulations of monopole magnetospheres of 
black holes which show the development of a powerful ultra- 
relativistic wind predominantly in the equatorial direction 
( Komissarov 2004b ). Being taken together, these two results 
suggest that divergence of magnetic field lines is a necessary 
condition for generating strong ultra-relativistic plasma out- 
flows from black holes. Moreover, given the fact that the 
Lorentz factor of the monopole wind reaches the value of 
W = 3 in the equatorial direction only at the distance of 
around r=20 (or 20r 9 , where r g = GM/c 2 , if one prefers this 
dimensional form) and much further away in the polar di- 
rection, one would not expected to find an MHD-driven col- 
limated ultra-relativistic outflow at distances smaller than 
several tens of r g . Here there may be present only rela- 
tivistic beams of particles accelerated by other mechanisms, 
e.g. electromagnetically. This makes the ongoing projects of 
studying the very bases (up to few tens of r g ) of astrophysi- 
cal jets via short- wavelength VLBI observations particularly 
interesting ijKrichbaum et aT~2 004 ) . 

The bottom right panel of fig0] shows the angular ve- 
locity of magnetic field lines, Q, normalised to the angular- 
velocity of the black hole ttn = <*/(*"+ + a 2 ), where r+ is the 
radius of the event horizon. In a steady state this parameter 
must be constant along magnetic field lines (e.g. Camenzind 
1986) which is what is seen in the plot. Thus, at t = 60 our 
solution is very close to a steady state at least for r < 20. 
Another interesting result seen in this plot is that £1 ~ 
0.5f2h, the value obtained in (Blandford & Znajc k~i977| for 
a slowly rotating black hole with monopole magnetic field. 
Finally, there is a sharp transition between the "rotating 
column" of magnetic field lines attached to the black hole 
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Figure 5. f2/f2/, (top panel) and Hj, (bottom panel) at r=5 and 
t=60. 

and the nonrotating "soup" of magnetic field lines which fail 
to enter the ergosphere (see also fig|HJ- This is most likely a 
discontinuity somewhat smeared by numerical diffusion and 
perhaps by the process of reconnection in the equatorial cur- 
rent sheet described above. 

The bottom left panel of fig^] shows the distribution of 
H<f,, the covariant azimuthal component of vector H intro- 
duced via 

Hi = *F U . 

In a steady-state 

V x H = J, 

where J is the electric current density (IKomissarov 2004al . 
Thus, H^, at the point (r, 9) gives us the total electric cur- 
rent flowing through the loop r =const originated from this 
point. In force-free steady-state solutions H,/, is also con- 
stant along magnetic field lines ( |Blandford fc Znajek 1977| 
IKomissarov 2004aTl . This is exactly what is seen in figure 0] 
thus confirming that at t = 60 the solution is almost force- 
free and very close to a steady state for r < 20. The distri- 
bution of H$ also exhibits a discontinuity at the boundary 
of the "rotating column" indicating a thin current sheet of 
return current. 

The split-monopole structure of magnetic field found 
in these simulations within the ergosphere is in conflict 
with the electrodynamic simulations where the field lines 
exhibit a sharp turning point in the equatorial plane 
(Komissarov 2004a). There seem to be only two possible 
reasons for this difference. First of all, neither the inertia 
nor the pressure of plasma particles are accounted for in the 
electrodynamic model. At first glance, this does not seem to 
be important as even in the current MHD simulations both 
are small compared to the mass-energy density and pressure 
of magnetic field almost everywhere. But not exactly every- 
where. In the equatorial current sheet the gas pressure dom- 
inates and plays a stabilizing role against otherwise quick 
reconnection of magnetic field lines. The second factor is 



the resistivity model. While the electrodynamic simulations 
(Komissarov 2004a I utilize the anisotropic resistivity based 
on the inverse Compton scattering of background photons, 
the only resistivity present in the MHD model is numerical. 
At this point we cannot exclude that explicit inclusion of 
physical resistivity will have a strong effect on the outcome 
of MHD simulations. This matter will have to be addressed 
in future studies. 



4 SUMMARY AND CONCLUSIONS 

We have carried out axisymmetric ideal MHD simulations 
of plasma flows in the case of a rotating black hole immersed 
into an initially uniform magnetic field described by the vac- 
uum solution of Wald {1974 1 . Our main intention was to 
verify the results of previous simulations (Koid e~et al. 20021 
IKoide 200311 claiming a rather significant role for the so- 
called MHD Penrose process of extracting the rotational en- 
ergy of black holes, at least in this particular case. Our sim- 
ulations show that this is true only for a short initial period 
during which a large fraction of magnetic field lines enter- 
ing the black hole ergosphere exhibit a turning point in the 
equatorial plane. Eventually, those field lines are pulled into 
the black hole and within the ergosphere the magnetic field 
acquires the split-monopole structure. After this transient 
phase the regions of negative hydrodynamic energy at infin- 
ity are no longer present in the ergosphere and the the MHD 
Penrose process ceases to operate altogether. 

The rotational energy yet continues to be ex- 
tracted via purely electromagnetic Blandford-Znajek pro- 
cess ( |Blandford fc Znajek 1977| IKomissarov 2004all . This 
energy extraction, however, is not followed by develop- 
ment of large-scale outflow from the black hole in sharp 
contrast with the results of previous MHD simulations 
for the monopole configuration of magnetic field lines 
( Komissarov 2004b ) . The only exception is a thin current 
sheet located at the interface between the "rotating col- 
umn" of magnetic field lines attached to the black hole and 
the "nonrotating soup" of field lines failing to enter the 
black hole ergosphere. The relatively slow outflow in this 
"sheath" is likely to be driven by gas pressure and is fed 
by hot plasma produced during reconnection events in the 
equatorial current sheet. Given the results of these and other 
recent MHD simulations (IKoide 20031 IKomissarov 2004bf it 
seems impossible to generate via the standard MHD mech- 
anism an outflow which becomes both ultrarelativistic and 
collimated already within the few first decades of the grav- 
itational radius from a black hole. Should relativistic as- 
trophysical jets continue to be seen at such small distances 
lIKrichbaum et al. 2004L one will have to look for other ex- 
planations of their origin. 

One of the main shortcoming of these simulations is 
the approximation of perfect conductivity. The only source 
of resistivity governing the process of magnetic reconnection 
in the developing current sheets is purely numerical. Future 
numerical models will have to include physical resistivity. 
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